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ABSTRACT 



Seasonal longshore flow patterns are examined at four points along 
the Oregon-Northern California coast. Summer and winter activity is 
examined as far seaward as 25 nautical miles and as deep as 200 meters. 
Long-term mean hydrographic data are used to determine geostrophic ve- 
locities. A nearshore baroclinic southward flow (^^20 cm/sec) is ob- 
served at each of the points during the summer. Winter currents are 
generally very small (<10 cm/sec) and largely barotropic in nature. 
Seasonal volume transports are presented; corrected velocity profiles 
are also presented based on data from moored current meters. Quali- 
tative explanations of the observed phenomena are considered. 
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I. INTRODUCTION 



A. DESCRIPTION OF AREA OF STUDY 

The area of the study, shown in Figure I, is the one studied by 
Bourke (1971). More specifically it is bounded by 41° and 46° North 
latitude. It extends from the coast of Oregon out to approximately 
25 nautical miles. Subsurface activity is examined down to a depth 
of 200 meters. 

The California Current flows through the study area and is the 
main influence on water movement. Significant variations from the 
California Current pattern may develop locally. These deviations 
may be attributed to winds, tides, the shallow topography or a com- 
bination of these elements. 

The topography is characteristic of the West Coast of the United 
States in that the shelf is very narrow and steep. Steepness general- 
ly increases southward as seen in Figure 2. Immediately north of 
Zone A is Astoria Canyon just off the mouth of the Columbia River. 

At the southern end of Zone B are located Heceta Bank and Stonewall 
Bank, features that shoal to depths less than 60 meters and may cause 
some deflection of the longshore current. The protrusion just south 
of Heceta Bank is Cape Blanco which divides Zone C from Zone D. 

A mixture of the Oyashio and the Kuroshio current systems supplies 
the California Current. The current is broad, slow, and shallow 
(--500 m) . It flows southward in a band roughly 300 miles wide at an 
average speed of 10 cm/sec (Dodimead, Favorite, and Hirano, 1963). A 
seasonal variation of the California Current is the Davidson Current. 
It develops off the Washington-Oregon coast in September, becomes well 
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Figure 1. The coastal boundary of the study area, 
(from Bourke , 1971) 
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Figure 2. Three-dimensional representation of the approximate topography of the study area. 



established by January and disappears by May. It develops inshore of 
the California Current and may reach 50 miles in breadth and attain 
speeds of 20 to 40 cm/sec (Schwartzlose , 1963). 

The Oregon coast is located approximately in the center of the 
zone of prevailing westerlies with local winds varying from northwest 
to southeast throughout most of the year. The barrier presented by 
the mountains of the Coast Range influences the general wind pattern, 
deflecting the winds so that they tend to align with the trend of the 
coast (Cooper, 1958). In regions where the mountains are low, deflec- 
tion is minimal and normal wind conditions exist. 

B. ORGANIZATION OF DATA 

The National Oceanographic Data Center (NODC) has on file all hy- 
drographic data taken in this region prior to 1966. Data from 1966 
to 1971 are those observations made by the Department of Oceanography 
of Oregon State University (OSU) . The NODC data are filed by 10^ Mars- 
den square numbers; number 157 encompasses the area of this study. 

Bourke employed a computer program which excluded all data taken 
more than 25 miles from shore. The data were then grouped by month so 
that summer and winter seasons could be determined. The summer season 
included all data sampled in July and August while the v;inter season 
included data from December through March. These data were found to be 
grouped about hydrographic lines frequently sampled by OSU. These lines 
originate off Astoria (AH line), Tillamook Head (TH line), Depoe Bay 
(DB line), Nevjport (NH line), Coos Bay (CH line), and off Brookings 
(BH line). Along hydrographic lines observations are made at regular 
intervals at 5, 15, 25, ... nautical miles offshore. Occasionally 
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additional stations are taken at three and ten miles. As a result, 
within the study area, most of the data were clustered about hydro- 
graphic lines oriented perpendicular to the coast and were further 
clustered at positions 5, 15, and 25 miles from the coast. 

With these points in mind Bourke established a grid system where- 
in average values of temperature and salinity could be computed for a 

\ 

specified area within the grid network. Average values were computed 
for depths of 0, 10, 20, 30, 50, 75, 100, 150, and 200 meters. 

The study area was divided into four zones as shown in Figure 3. 
Factors affecting the temperature and salinity were assumed to be uni- 
form or consistent within each zone. The zone boundaries were con- 
structed after having considered the shape of the coastline, location 
of coastal rivers, bottom topography, and distribution of data. Zone 
A extends from 46^ 16* N to 45° 38* N, Zone B from 45° 38* N to 44° 19* 
N, Zone C from 44° 19* N to 42°47*N, and Zone D from 42° 47 * N to 41° 

41* N. Each zone encompasses at least one hydrographic line. 

Although the average values apply to a given area, they were as- 
signed to a specific point for calculation purposes. These artificial 
stations provided a network for the computation of geostrophic veloc- 
ity and transport within the study area. 

C. PREVIOUS STUDIES 

Bourke analyzed the data used in this thesis as regards temperature 
and salinity structure. In general, he found that the surface waters 
were warmer in summer than in winter. Summer temperatures were 2 to 5°C 
warmer than winter temperatures with the largest increases occuring off 
Northern Oregon. Surface temperatures increased with distance offshore 
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Figure 3. Schematic diagram of data grid showing zone boundaries. 

Hydrographic lines are shown as dotted lines. (from 
Bourke, 1971) 
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in both seasons. The offshore temperature gradient was 4 times larger 
in summer than winter. In the longshore direction surface temperature 
increased toward the north in summer while the reverse was true in win- 
ter. The temperature gradient was nearly the same in both seasons. 

Surface waters were 0.5 to 1.0 ppt more saline in summer than in 
winter. Surface salinities decreased with distance from shore in sum- 
mer. In winter surface salinities increased with distance from shore. 
Salinity gradients in the x (offshore) direction were identical in sum- 
mer and winter. In both seasons surface salinities increased toward 
the south. The longshore gradient, however, was twice as large in 
summer as in winter. 

Variability in temperature, salinity, and density with distance 
from shore was significant in summer but relatively unimportant in 
winter. A strong thermocline was present in summer from the surface to 
30 meters. Winter surface waters were isothermal to 50 meters. A 
strong halocline was observed in both seasons extending to .75 meters 
in the summer and 100 meters in winter. In both seasons at 200 meters 
the temperature and salinity were essentially constant throughout the 
study area. At 200 meters sun^ner temperatures and salinities were 
cooler and more saline than winter values due to the strong coastal 
upwelling occuring during the summer months. Summer values were 
6.8 ± 0.2°C and 33.98 ± 0.01 ppt; in winter these wore 7.7 db 0.1°C 
and 33.90 ± 0.5 ppt. 

The frontal layer (permanent pycnocline) , defined by Collins 
(1964) to be bounded by the 25.5 and 26.0 surfaces, rises sharply 

shoreward in summer, but is essentially flat in winter. Its depth in 
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winter is about 100 meters. During the summer it slopes from a depth 
of about 50 meters at 25 miles offshore to a depth of about 20 meters 
at 3 miles offshore. In Zone D it breaks the surface at about 7 miles 
offshore . 

Pavlova (1966) performed a study of the seasonal variations in the 
California Current and Countercurrent in latitudes 20^ to 40^N. Using 
hydrographic data averaged over 35 years she computed geostrophic ve- 
locities and transports relative to 1500 decibars. She also mapped 
the dynamic topography of the California Current. 

In developing a characteristic profile of the California Current 
and Countercurrent, she chose a cross section oriented perpendicular 
to the coast around 32^ north latitude. A northward countercurrent 
is observed at depth throughout the year. Except during the spring, 
a surface countercurrent is also observed throughout the year. Max- 
imum countercurrent velocities (8 to 9 cm/sec) are reached in the sum- 
mer and autumn at a depth that varies in the course of the year: 
from 200 to 300 meters at the beginning of the summer and in the au- 
tumn, to 100 meters at midsummer. On an annual average, however, 
the countercurrent at depth is practically a single and uniform flow 
from 25^ N to 35^ N and farther. 

The strength of the countercurrent is related to the California 
Current. Maximum velocities for both currents are reached in summer 
and in autumn when the northerly geostrophic wind is strongest and up- 
welling is most intense along the coast. 

Pavlova also notes that in the upper layers there is an upward 
and in the lower layers a downward isopycnic slope from the ocean 
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towards the coast, so that at intermediate depths (100 to 250 meters) 
a zone of divergence of the isopycnals forms in relation to the posi- 
tion of the profile and the time of year. Her data indicate that the 
countercurrent reaches maximum development in this zone of divergence, 
which in August is centered on the 26.0 isopycnic surface. 
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II. THEORY 



A. INTRODUCTION 

It is assumed that the long-term mean current is composed of two 
parts: a baroclinic component arising from the distribution of the in- 
ternal mass field and a barotropic component arising from the slope 
of the sea surface. Geostrophy is assumed in order to determine the 
vertical velocity profile and this is transformed to an absolute pro- 
file by the use of data from moored current meters, 

B. GEOSTROPHIC VELOCITY 

The theoretical development is presented in the left-handed co- 
ordinate system where the z axis is positive downward and the y axis 
is positive northward. The components of velocity, u and v, are di- 
rected along the x and y axes, respectively. Pressure is denoted by 
p, density by p, and the eddy viscosity coefficient by A. Other de- 
finitions to be used are specific volume, o', defined as the recipro- 
cal of density and the Coriolis parameter, f, defined as 2Qsin0, where 
0 is latitude and Q is the angular velocity of the earth. 

Since north-south (meridional) velocities are being examined, 
the appropriate equation of motion describes the east-west (zonal) 
acceleration : 



This equation states that the frictional forces present can be re- 




The frictional term is chosen to have the form: 




( 1 ) 
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presented as shearing stresses. Furthermore it is assumed that the 
lateral shearing stresses are minor in comparison with the vertical 
shearing stresses. The equation of motion can then be written as 



where the V notation has been dropped on the vertical eddy viscosity 
coefficient . 

This equation is simplified by two assumptions based on the fol- 
lowing considerations : 

(a) The data consist of long-term monthly averages of salinity 
and temperature. These values are further averaged over a 
seasonal period of tw^o months or more, so that a smoothed 
seasonal description of the study area may be presented. 

(b) Frictional effects, although present, are confined to rela- 
tively thin boundary layers (of the order of 20 meters) at 
the sea surface and along the bottom of the shelf and slope 
regimes. These effects exert an insignificant influence on 
the geostrophic interior of the area. 

In support of (b) the following calculations are presented for 
Zone B where sufficient current measurements were available. The 
average wind speed for July and August is 6.0 m/sec blowing roughly 
from the north most of these two months (Fleet Numerical Weather Cen 
ter, Monterey, California). Based on a v/ind speed of 6.0 m/sec, Neu 
man’s (1966) value for the "effective” eddy viscosity coefficient, A 
is 160 gm/cm-sec. The following current measurements were collected 
by OSU seven nautical miles offshore of Depoe Bay, Oregon, in 100 



_ du 

dt 



fv - q|^ + q|-(A— ) (2) 

OX oz 7)z 



meters of water, from 23 June to 10 July 1968: 
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Depth of sensor 

(e.} 



u component* 
(cm/ sec) 



25 


0.9 


50 


2.1 


75 


4.8 



* Velocities are averages based on common record 
lengths, (adapted from Pillsbury, et al., 1970) 

If A is assumed constant and the above data are used to evaluate 
3^u/5z^ numerically, one obtains a value of 2.4 x 10 (cm-sec) . 

The frictional wind stress term is, therefore, of the order of 10 
while the Coriolis term, fv, is of the order of 10 . 

By reason of (a) above, du/dt = 0. By reason of (b) , 

OL A 3^u/3z^ 0, so that the equation for the x component of motion 

may be expressed as 

fv = O' r-^ 

3x. (3) 

It is recognized that the time averaging process which eliminates 
the time dependence in equation (2) also reinforces the stress terms 
of equation (1). This result is deemed negligible and equation (3) 
is used as the governing equation for the distribution of the meri- 
dional velocity field. 

For small pressure differences, Ap , and small distances, Az and 
Ax, one may express the vertical and horizontal pressure gradients as 

Ap = pg Az 

and Ap = (Sp/3x) Ax. 

Eliminating Ap yields 

g Az = (3p/3x) Ax. 

Rewriting the equation leaves 
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Ax 



a 



Sx. 



In the limit that and Az approach zero, while the pressure remains 
constant, one obtains 



“ 9 X. (A) 

The pressure gradient in equation (3) can now be replaced by the zonal 
slope of an isobaric surface. The concept of a slope, however, implies 
the existence of an absolute reference surface (i.e. a level surface), 
a surface frequently unknown in oceanographic work. Thus the idea of 
a relative slope is introduced. Given the vertical distances, h and 
h , between isobaric surfaces at stations A and B, respectively, the 

D 

slope of surface p^ relative to a deeper surface P 3 is given by 

(|^) = ~ \ 



where L is the distance between stations and station A is east of sta- 
tion B as illustrated in Figure 4. 




Figure 4. Determination of the slope of the pi 
surface relative to the p^ surface. 
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The evaluation of and is accomplished utilizing the concept 
of geopotential as described by Sverdrup, Fleming and Johnson (1942). 
The result of their discussion states that a geometrical distance, z, 
is proportional to the pressure integral of the specific volume of the 
water column, 

z = “ J O' dp. 

A convenient, but arbitrary, choice of units requires that g be spe- 
cified in m/sec^, p in decibars (10^ dynes/cm^) and O' in cm^ /gm. The 
depth, z, is then specified in meters. Applying this technique to 
equation (5) yields 






( 6 ) 



..P2 ,.P2 

J ' 

Pi Pi 

Substituting equations (3) and (4) into equation (6) gives the rela- 
tive current perpendicular to the plane between stations A and B, 



Vi - Vs = 



10 

fL 



r 1 

(J O' dp)j^ - (j O' dp)J 
Pi Pi 

The integral over pressure is the geopotential distance, D, expressed 
in dynamic meters. According to Sverdrup, et al., D is defined as the 
sum of Dq and AD, where 



(7) 



Do = Ds - D, = J O' dp 

Pi 

and A D = J 6 dp. 

Pi 

The standard geopotential distance, Dq , results from integrating over 
a "standard" water column of salinity equal to 35 ppt and temperature 
equal to 0°C. The geopotential anomaly, 4D , results from integrating 
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over a water column whose salinity and temperature structure are de- 
fined solely in terras of deviations from the structure of the ’’standard" 
water column. The symbol, 6, denotes the specific volume anomaly. Re- 
writing equation (7) gives 

Vi - Vg = ~[(J 6 dp)^ - (J^'s dp)J. (8) 

Pi Pi 

Since the ocean is very nearly in hydrostatic equilibrium and 
since the decibar was chosen as the dimension of pressure, dp may be 
replaced by dz , where z is measured in meters. As a working approach, 
a level sea surface is assumed so that v^ vanishes. With this assump- 
tion, V 2 is necessarily a baroclinic velocity only. It is caused 
strictly by the distribution of the internal mass field. There is no 
barotropic component since there is no slope of the sea surface. The 
integration is performed from the sea surface to the depth of interest, 
therefore 



v(z) = V 2 = 



10 

fL 



z z 

[(J 6 - (J 6 dz)g_^ 



(9) 



With this formula v(z) is rendered in meters per second. By speci- 
fying L in kilometers and omitting the factor of 10 in the numerator, 
v(z) may be calculated in centimeters per second. An observer, look- 
ing at the plane of the stations with station A to the right of sta- 
tion B, will note that a positive current flows away from him while a 
negative current flows toward him. 

By means of a fairly uncomplicated, but lengthy, computer program, 
the specific volume anomaly may be determined by the appropriate poly- 
nomial (Knudsen, 1.901). Numerical integration to successive depths 
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generates the baroclinic velocity profile. However, the profile is 
still relative to the sea surface which has been assumed level. As 
it turns out this is usually a poor reference level and the calculated 
shear must somehow be fixed at a point to provide a more realistic dis- 
tribution of velocity. Traditionally this has been done by assuming 
that horizontal motion ceases at some great depth, usually around 
1500 decibars (Defant, 1961). The profile is then adjusted to show a 
zero current at 1500 db. In shallower areas a more positive approach 
is to measure the current at some depth (or depths) and then adjust 
the baroclinic profile to make the best fit to the data. In effect, 
this amounts to adding a barotropic component to the velocity profile. 
Tliis latter technique is equivalent to a physical determination of the 
reference level. Both techniques allow the computation of the sea 
surface slope once the shear has been altered. In the computer pro- 
gram, the user specifies a reference level, and the vertical velocity 
profile is adjusted so that the current is zero at this depth. Know- 
ledge of the correct reference level is imperative if the profiles 
are to provide useful information. 

C . VOLUME TRANSPORT 

When describing oceanic circulation, a quantity of interest is 
the volume transport between two stations. Consider a plane section 
bounded on either side by stations A and B, at the top by the sea sur- 
face and at the bottom by some isobath,, d. The volume transport is a 
measure of the amount of water passing through this plane section per 
unit time. Given the velocity as a function of depth, the volume 
transport is 
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d 

T = L J V dz, (10) 
o 

Using the computer program mentioned previously, it is a simple matter 
to insert the adjusted velocity values into this equation and perform 
the indicated integration. Equation (10) is preferred to the often 
quoted equation, 



T = J (ADg - AD^) dz. (11) 
o 

This equation is developed by making the following substitution for v 
in equation (10) : 



12 5(ZM)) 10 /SD - AD 

f Bx ^ f S 



where D and D are the geopotential anomalies at stations A and B, 

A 15 

respectively, in the configuration described earlier. Application of 
equation (11) requires that the velocity must vanish at some depth. 

No such requirement is imposed on equation (10). Using the adjusted 
profile, transport in water layers can be calculated and then summed 
to give total volume transport. To a sufficient degree of accuracy, the 
mass transport may be determined by multiplying the volume transport 
by an average density value. 



D. THERMAL WIND RELATIONSHIP 

The thermal wind equation, sometimes written as 

g=-fEs (12) 

will be referred to in the discussion of procedure. The simplifying 
conditions needed for the development of this expression are given by 
Pillsbury (1972) in a ver>" straight-f on\/'ard derivation. Using the 
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same coordinate system as before and with only the pressure and Cor- 
iolis forces acting in the east-west direction, the equation of motion 
becomes 



p f V 



3x 



Using the hydrostatic equation, 



P g = 



- ^ 



Sz 



the pressure term is eliminated by cross differentiation to leave 

5v _ __g ^ V ^ 

Sz fp Bx p 3z . (13) 

In the study area, v is of the order of 10 cm/sec making 

^ ^ ^ 10 
fp p 



The slope of the isopycnal. 



.5zs _ ^ 5o /5x 
ox p 3p/dz 



is not often less than 10 ^ in the study area (Pillsbury, 1972). 

Therefore the ratio. 



_£ ^ / z le. 

fp 3x p 5z 

is of the order of 10 ^ . Thus the second term on the right-hand side 
of (13) may be neglected, and the vertical shear may be written as 

5v _ g 5P 
Sz fp Sx . (14) 

Now the static stability is given by 



E = i|£ 

p oz 



and the product of the isopycnal slope and the static stability is 
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given by 



E s = 

p OX . 

Substituting this expression in (14) leaves equation (12), 



5v 

3z 




s. 
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Ill, PROCEDURE 



A. INTRODUCTION 

The original objective of this study was to compute relative sea- 
sonal velocity profiles in the four zones and then adjust these pro- 
files to agree with observed current measurements. It was thought that 
this procedure would develop a reliable long-term mean picture of the 
water circulation over the Oregon and Northern California shelf-slope 
area. With the corrected velocity profiles, meaningful longshore trans- 
port computations could be made. 

These plans had to be modified due to the lack of sufficient cur- 
rent observations. The best measurements available were collected by 
Oregon State University* s Department of Oceanography. These were 
taken in coastal waters off Depoe Bay and Newport, Oregon, in 1966, 1967, 
1968, and 1969. They were used to adjust the velocity profiles of 
Zone B. Other measurements were taken by the University of Washington *s 
Department of Oceanography. It was hoped that the latter could be 
used in the profiles of Zone A since they were collected near the mouth 
of the Columbia River. This was not the case, hov^?ever, since the mea- 
surements that were collected during the appropriate months, were col- 
lected at depths deeper than the available hydrographic data. No 
current data were located for Zones C and D. In addition scarcity 
of winter hydrographic data in Zone D precluded the formation of sta- 
tistically significant average values for that season. 

Consequently the results presented are not absolute, but do illus- 
trate a fairly accurate description of the variations in relative water 
movements from summer to winter and alongshore from north to south. 
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B. METHOD 



A computer program was obtained to compute the geostrophic veloci- 
ties via equation (9) as developed in section II. The use of geostro- 
phy is based on the assumption that the frictional boundary layers at 
the bottom and at the sea surface are small in comparison with the depth 
of most of the area. The wind, which is an important consideration in 
this region, is assumed to control only the barotropic component of 
the flow due to its influence on the sea surface slope. 

The problem of wind stress has been reasonably refuted by the nu- 
merical example cited in the theoretical development, but some assump- 
tions were involved there that need to be examined. The 1968 current 
data used in the determination of S^u/9z^ were collected in a period 
that did not fall entirely within the months of July and August, the 
period used for geostrophic calculations. It was assumed that during 
these months the currents did not vary significantly from the observed 
currents during the time 23 June to 10 July 1968. Also specific wind 
statistics were not available for the station at the time the current 
measurements were being recorded. Consequently, the wind speed at 
the time of observation was assumed to be the long-term mean wind 
speed as tabulated by FM-JC for the months of July and August (July 
mean, 12.6 kt.; August mean, 11.0 kt.). 

The wind speed was the determining factor in choosing a value for 
the eddy coefficient, A. In order to increase the magnitude of the 
frictional term to the point where it would be significant with re- 
spect to the Coriolis term, it is estimated that the mean wind speed 
from 23 June to 10. July 1968 would have to have been about 20 knots. 
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The applicable eddy coefficient would then have a value of about 580 
gm/cm-sec instead of the value 160 gm/cm-sec. The frictional terra 
would therefore amount to about 10 percent of the Coriolis terra and 
probably could not be neglected in treating the equation of motion. 

The assumption that mean wind conditions prevailed at the time 
of data collection might be open to question. Nevertheless, one wishes 
to have confidence in long-term mean data when applying it to any 
one particular time segment. If the number of observations of wind 
speed is large and these observations are concentrated about the mean, 
then the assumption of mean conditions is probably adequate. This 
is the case in the study area during the summer. 

The irregular topography of the area will complicate the geostro- 
phic flow pattern but will not affect the condition of geostrophy 
per se except in the very shallow regimes. Dynamically the topo- 
graphic anomalies do not appear to influence the hydrographic data. 

For a given depth, temperature and salinity do not vary significantly 
in the longshore direction within the zone boundaries. This fact is 
one of the points considered by Bourke when he first established the 
grid system of this area. 
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IV. RESULTS 



A. VELOCITY FIELDS IN SUMMER 

The most pronounced results are presented in Figures 5 through 8. 
These profiles represent geostrophic meridional velocities that were 
computed relative to the deepest common data point between stations. 
They do not represent accurate absolute values, but they are reliable 
in terms of the distribution of vertical shear. Examined in this 
light, the profiles illustrate some interesting patterns. 

The diagram for Zone A during the summer shows a surface jet start- 
ing at 13 to 15 nautical miles offshore and extending to the seaward 
limit of the data (25 n. mi.). This jet is the result of the dis- 
charge from the Columbia River. Using theory developed for surface 
discharges from ocean outfalls, Bourke has shown that the temperature- 
salinity characteristics of the core of the Columbia River plume remain 
intact to a distance of 10 to 15 miles from the mouth of the river. 

The jet in Zone A profile, therefore, represents a southward flowing 
mass of warm and relatively fresh water that has been discharged 
from the jetties at the mouth of the Columbia. The averaged hydro- 
graphic data support this analysis since cold saline oceanic water is 
found on both sides of the jet which is characterized by its warm, 
relatively fresh flow. 

The Columbia River discharge exerts a minimal influence on Zones 
B and C. In these zones the distinctive temperature-salinity of the 
core of the Columbia River plume is not present and oceanic water pre- 
dominates. The Zone B diagram also shows a surface jet located in the 
interior of the zone. Other barotropic data could shift the location 
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Figure 5. Average summer meridional velocity 
field for Zone A in cm/sec. 




Figure 6. Average summer meridional velocity 
field for Zone B in cm/sec. 
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Figure 7. Average summer meridional velocity 
field for Zone C in cm/sec. 




Figure 8. Average summer meridional velocity 
field for Zone D in cm/sec. 
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of the jet east or west, but 10 nautical miles is a reasonable posi- 
tion. This jet flow is a characteristic predicted by a numerical model 
of this area developed by 0*Brien and Hurlburt (1972). They used 
a model consisting of two Viomogenous layers coupled together by shear- 
ing stresses but not by mixing. The system is driven by a meridional 
wind stress. Their results predict a southerly current centered 10 
kilometers offshore, in reasonable agreement with the value of 10 
nautical miles (20 km.) predicted here. 

The Zone C profile is based on less data than Zones A and B. Of 
primary concern here is the absence of sufficient hydrographic data 
to construct long-term mean values at a position less than 5 miles 
off the coast. As a result any isotachs drawn shoreward of 10 miles 
are speculation. Nevertheless, the profile indiccites that the jet, 
if it is the same one from Zone B, has moved seaward. It is note- 
worthy that here, as in the other zones, the jet remains positioned 
above the steeper portion of the permanent pycnocline as described by 
Bourke. 

What might be happening in Zone C is that the northerly wind 
along the coast is being constricted and accelerated as it passes Cape 
Blanco. The physical result is an area of uniformly high winds and 
decreased pressure. Consequently one observes more intense upwelling 
in this area than in Zone B. 

Zone D, located just south of Cape Blanco, appears to exhibit 
the effect dramatically, Upwelling has intensified as shown by the 
presence of cold saline water farther out to sea and by the permanent 
pycnocline breaking the surface at about 7 miles offshore as mentioned 



32 



earlier. 



As discussed beforehand these profiles were computed relative 
to the deepest common data point between stations. One cannot expect 
this method to be very accurate in an absolute sense, but it is more 
realistic than assuming surface currents are zero. As a check on 
the procedure thus far, Figures 9 and 10 are presented. These are pro- 
files for Zone B that have been corrected using current measurements 
from moored current meters. The data, collected and processed by OSU, 
are listed in Tables I and II. 

The measurements from 1966 indicate a vertical shear of 4.7 x 10 
sec between the depths of 20 and 60 meters at 5 nautical miles 
offshore. Using mean hydrographic values in the thermal wind rela- 
tion yields a vertical shear of about 4.2 x 10 sec between sta- 
tions B-05 and B-15 over an D.ntexrval of depth centered at 40 meters. 

The 1969 observations show a vertical shear of 2,6 x 10 " sec be- 

tween the depths of 20 and 80 meters at 15 nautical miles offshore. 
Between stations B-15 and B-25, the average hydrographic data in the 
thermal wind relation gives a vertical shear of 1.6 x 10 sec 
over an interval of depth centered at 50 meters. Tlie agreement be- 
tween the measured shear and the predicted geos trophic shear is re- 
assuring and lends a little more credence to the adjusted profiles. 

Consistent with the profiles discussed previously, the adjusted 
profiles (Figures 9 and 10) confirm the existence of the surface jet. 
The presence of a subsurface countercurrent is indicated by the 1966 
profile. These results have also been noted by Mooers , Collins, and 
Smith (1972). Their data included that of Table I and a line of 
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Figure 9. Summer meridional velocity field for Zone B 
adjusted with 1966 current observations. 




Figure 10. Summer meridional velocity field for Zone B 
adjusted with 1969 current observations. 



34 



DEPTH (nn] _ DEPTH [m) 



Table I.* Mean Current Statistics, 1966 

(based on hourly averages for common 
record length of 339 hours) 



Depoe 

Bay 

Station 


Distance 
offshore 
(n mi) 


Water 

depth 

(m) 


Sensor 

depth 

(m) 


Eastward 
c omp one n t ** 
(cm/ sec) 


Northward 
component*' 
(cm/ sec) 


5 


5 


80 


20 


-2.2 i 11 


-17.9 ± 12 


5 


5 


80 


60 


2.8 db 6 


0.7 ± 9 


10 


10 


140 


20 


4.5 ± 10 


-12.1 ± 9 


15 


15 


200 


60 


6.0 ± 6 


7.9 ± 13 



Adapted fromMooers, Collins, and Smith (1972). 



•kic Velocities are presented as the mean value ± the standard 
deviation . 



Table II.- Mean Current Statistics, 1969 

(based on total record lengths 
acquired in August and September) 



Station 
(Newport or 
Depoe Bay) 


Distance 
offshore 
(n mi) 


V7ater 

depth 

(m) 


Sensor 

depth 

(m) 


Eastward 
component** 
(cm/ sec) 


Northward 
component-' 
(cm/ sec) 


NH-3 


3 


50 


20 


0.7 ± 7 


*4 . 4 i 16 


DB-7 


7 


100 


40 


-6.1 ± 11 


-21.1 ± 19 


NH-15 


15 


100 


20 


-9.9 ± 16 


-18.1 16 


NH-15 


15 


100 


80 


-3.5 ± 14 


-2.8 ± 12 


* Adapted 


from Huyer, 


et al. 


(1971). 







-- Velocities are presented as the mean value ± the standard 
deviation. 
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hydrographic stations taken in late July, 1966, off Newport, Oregon. 
The velocity field presented by Mooers, et al. differs only slightly 
from the profile presented here based on long-term mean hydrographic 
data and the data from Table I. The 1966 profile centers the south- 
ward flowing jet at about 5 nautical miles (10 km.), in better 
agreement with the results of 0*Brien and Hurlburt. The northward 
undercurrent apparently does not break the surface as suggested in 
Figure 9. Mooers, et al. indicate that in 1966 the surface flow 
was all southerly at least as far as 60 miles offshore. They report 
that the undercurrent may have come within about 10 meters of the sur 
face that summer. 

B. VELOCITY FIELDS IN WINTER 

The winter profiles are presented in Figures 11 through 13. 

Prior to assuming a reference level at depth, the relative (baroclin- 
ic) current values were observed to be very small, illustrating the 
importance of the barotropic component during winter. A zonal sea 
surface slope of one millimeter per kilometer, corresponding to a cur 
rent speed of the order of 10 cm/sec, would be sufficient to reverse 
the flow at any point in any of the zones. 

In spite of such small geostrophic velocities a description of 
a characteristic flow pattern is attempted. It is felt that the long 
term nature of the data used in this analysis provides a firmer 
ground for such an attempt. 

Nearshore surface waters in Zones A and B are responding to the 
wind stress while the offshore and deeper waters are moving under the 
influence of the large scale oceanic circulation. Current velocities 
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Figure 11. Average winter meridional velocity 
field for Zone A in cm/sec. 




Figure 12. Average winter meridional velocity 
field for Zone B in cm/sec. 
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Figure 13. Average winter meridional velocity 
field for Zone C in cm/sec. 
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are greatly reduced from the summer values due to the nearly level 
isopycnal surfaces. Runoff and water piled up onshore by the x^7ind 
contribute a northerly barotropic component to the circulation. 

The University of Washington has made a few current observations 
during the winter season at moorings near the Astoria Canyon (Hopkins, 
1971). Unfortunately, these observations were not sufficiently dis- 
tributed spatially to allow an absolute rendition of the velocity 
field. In each of these cases, however, the mean northerly component 
was only slightly larger (in the sense of being more positive) than 
the predicted geostrophic component. This fact illustrates tx^/o points: 
that the fresh water runoff and the water piled up onshore by the pre- 
vailing x^/ind govern the barotropic component of the velocity and that 
this barotropic component is not very large. 

Cc VOLUME TRANSPORT 

Volume transports x^ere computed using equation (10). As in the 
velocity computations, the deepest depth of the data common to the 
two stations x^as used as the reference level. Again it is cautioned 
that this is not necessarily the correct level of no horizontal motion, 
but it is better than a level sea surface approximation. 

The results of this endeavor are presented in Tables III and IV. 

As expected, the transports are uniformly small, the winter transports 
being roughly an order of magnitude smaller than the summer transports. 
(This observation is probably correct; hox^ever, it is not completely 
justified from an inspection of the tables alone. The winter data are 
generally shallox^7er than the summer data, thereby eroding a basis for 
comparison.) Reflecting the structure of the velocity fields, the trans- 
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Table 111. Geostrophic Meridional Transports (Summer) 



Stations 


Volume Transport 
(Sverdrup) 


Direction 


Reference Level 
(m) 


A-03/A-05 


0.010 


S 


50 


A-05/A-15 


0.067 


S 


75 


A-15/A-25 


0.411 


S 


200 


B-03/B-05 


0.018 


S 


75 


B-05/B-15 


0.201 


S 


150 


B-15/B-25 


0.193 


S 


200 


C-05/C-15 


0.119 


S 


75 


C-15/C-25 


0.287 


S 


200 


D-05/D-15 


0.134 


S 


100 


D-15/D-25 


0.065 


S 


200 



Table IV 


• Geostrophic Meridional Transports 


inter) 


Stations 


Volume Transport 
(Sverdrup) 


Direction Reference Level 

(m) 


A-03/A-05 


0.001 


N 


50 


A-05/A-15 


0.006 


S 


75 


A-15/A-25 


0.198 


S 


200 


B-03/B--05 


0.001 


N 


50 


B-05/B-15 


0.002 


S 


75 


B-15/B-25 


0.023 


S 


75 


C-05/C-15 


0.034 


s 


75 


C-15/C-25 


0.043 


s 


200 
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port appears to be southward in both summer and winter. While this 
result is probably true for the summer months, it is uncertain that 
this is the case during the winter months. Deeper data and an accu- 
rate determination of a reference level could drastically change the 
winter transports. The small nearshore northward transport observed 
in Zones A and B during the winter shows the effect of the wind in 
moving the surface waters. 
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V. CONCLUSIONS 



The fact that long-term data were used in this discussion allows 
more general statements to be made than could be made othen^7ise• 

The surface jet, for example, is shown to be a regular phenomenon of 
the summer season. Its position above the steepest portion of the up- 
turned pycnocline and its absence in winter tie it to the occurrence 
of coastal upwe.lling. The fact that it is a dominant feature of the 
average summer velocity fields indicates that it must be a steady 
state event or at least an average of a series of consistent events. 
Furthermore, the jet is shown to be a characteristic of most of the 
Oregon and Northern California coastline, 

Mooers, et al, have shown that a two-layered baroclinic model 
of the flow over the Oregon shelf and slope is at least partially 
valid. They postulate that the barotropic and baroclinic radii of 
deformation should provide fundamental length scales. The barotro- 
pic radius of deformation is defined as the distance a shallow water 

\h —1 

wave would travel in an inertial period: R, . = (gh) f 

^ barotropic 

The baroclinic radius of deformation is defined as the distance an 

internal wave would travel along the interface between a shallow 

upper layer and a deep lower layer of water in the same inertial 

period: R, . = (gD Ap/p)^ f~” , The estimated values arc 

^ baroclinic 

R, . = 2000 kilometers and R, -i . . =20 kilometers. By 

barotropic baroclinic 

plotting the depth of the pycnocline as a function of offshore dis- 
tance on semi-logarithmic axes, Mooers, et al. confirmed the expo- 
nential character of the *’interf ace” and determined its radius of 
deformation to be about 40 kilometers over the shelf and 70 kilometers 
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over the slope. By comparison then, the water is deemed sufficiently 
stratified to warrant use of a two-layered (baroclinic) model, al- 
though different length scales are in effect over the shelf and the 
slope regimes. 

Without becoming involved in cause and effect arguments, it is 
asserted that the longshore northerly wind, acting on a relatively 
thin surface layer, moves water offshore. This lighter surface 
water moves seaward without significantly disturbing the lower hea- 
vier water and causes a pressure gradient force directed onshore, 
which in turn gives rise to a southv/ard flow. Because the southerly 
flow diminishes seaward, vorticity is decreased (f = constant, 

Sv/Sx < 0) , If potential vorticity is to be conserved (Stommel, 1966), 
then the thickness of the surface layer must decrease in the vicinity 
of the horizontal shear. Thus one observes the surface jet, the 
sloping pycnocline, and upwelling as interrelated phenomena. 

Another feature that appears to be related to the upwelling pro- 
cess is the poleward flowing undercurrent. The presence of this sub- 
surface flow has been reported by Hooers, et al. and by other workers, 
but no theory has been developed to explain it. It is undeniably 
present in the 1966 profile of Zone B (Figure 9) and, given more 
shallow reference levels for geostrophic computations, it is not dif- 
ficult to imagine its presence at the seaward edge of the profiles 
presented in Figures 5 through 8. 

It is reasonable to expect that the inclined isopycnals nearshore 
are capable of producing an offshore pressure gradient force if the 
sea surface does not develop enough slope to counteract it. Under 
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the proper conditions therefore, one may observe the resulting north- 
ward baroclinic flow located beneath the pycnocline and offshore. In 
the shallow shelf regime, the effect is masked by shearing stresses 
with the bottom and with the nearshore surface jet. Tliis may explain 
why the subsurface flow is found during some summers and not found 
during others. 

These conclusions are compatible with those of Mooers, et al. They 
are in basic agreement with the results of O’Brien and Hurlburt, but 
here there is a notable discrepancy. Upwelling is kno\-m to consist of 
a transient and a steady state condition. This discussion cannot pro- 
vide an insight into any transient solutions, but it does indicate 
that the nearshore jet and possibly the subsurface flow develop into 
steady state features. O’Brien and Hurlburt, in obtaining the near- 
shore jet, do not .permit this result. Indeed their model does not 
have a steady state solution so that a discussion of such fundamental 
differences may not be valid. 
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VI. SUMMARY 



The features observed in the study area are attributable to gen- 
eral oceanic circulation and to local wind and climate. The California 
Current is the main influence on water movement over the shelf and 
slope. In general, the flow is southv/ard regardless of the season. 
Summer winds are northwesterly while winter winds are southeasterly. 
These winds govern the sea surface slope and thereby determine the 
barotropic flov; during the different seasons. Surface waters re- 
spond to the wind field and, in general, flow parallel to it. Cli- 
mate along the Oregon coast provides plentiful fresh water runoff 
during the winter and considerably less runoff during the summer. 

The influence of the Columbia River discharge that is felt so strongly 
in Zone A during the summer is much less evident during the winter 
due mainly to the shift in wind. 

As a result of these characteristics more activity is observed 

in the summer than in the winter. The summer wind pattern stimulates 

the process of upwelling which causes tilting of the isopycnals which 

in turn generates the baroclinic flow. The winter wind pattern piles 

\ 

up water onshore but apparently to no great extent. Consequently 
there is little barotropic flow and even less baroclinic flow. 

Summer circulation is characterized by the upwelling of cold 
saline water next to the coastline. Typically the surface salinity 
increases and the surface temperature decreases shoreward and south- 
ward indicating more intense upwelling in the lower latitudes of the 
area. Also observed during the summer is a nearshore surface jet 
flowing southward. The surface jet occurs in the region above the 
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steepest portion of the permanent pycnocline where the near surface 
isopycnals have relatively small slope. In Zone A the jet appears 
to be caused by the discharge of the Columbia River, whereas in the 
other three zones it appears to be intimately related to the upwelling 
process. A subsurface poleward flow is observed but is not conclu* 
sively documented as a regular occurrence. This flow develops below 
the pycnocline in a region where the isopycnals are tilted sharply 
under the influence of upwelling. 

Winter flow is characterized by small barotropic and baroclinic 
components. Surface salinities and temperatures increase seaward 
and southward indicating the greater effect of runoff in the higher 
latitudes of the study area. The very small velocities presented 

in Figures 11 through 13 suggest that no distinctive longsliore flow 
pattern exists during the winter season, altliough the surface waters 
in Zones A and B appear to be responding to the wind field. 

In both summer and winter the volume transport is mostly south- 
ward. Winter transports are consistently less than summer transports 
and they are usually less by an order of magnitude. 
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